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Abstract 



The ground-state properties of the recent reported proton emitter Tm have been studied 
within the axially or triaxially deformed relativistic mean field (RMF) approaches, in which the 
pairing correlation is taken into account by the BCS-method with a constant pairing gap. It is 
found that triaxiality and pairing correlations play important roles in reproducing the experimental 
one proton separation energy. The single-particle level, the proton emission orbit, the deformation 
parameters f3 = 0.22 and 7 = 28.98° and the corresponding spectroscopic factor for 145 Tm in the 
triaxial RMF calculation are given as well. 

PACS numbers: 21.10.-k, 21.10.Dr, 21.10.Jx, 21.60.-n 
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I. INTRODUCTION 



Proton-rich nuclei display many interesting structural properties which are important 
both for nuclear physics and astrophysics. These nuclei are characterized by exotic decay 
modes, such as the direct emission of charged particles from the ground state, and /3-decays 
with large Q-values. The decay via direct proton emission provides a unique insight into the 
structure of nuclei beyond the drip line limit. The evolution of the single-particle structure, 
nuclear shapes and masses can be deduced from measured properties of proton emission [l| . 

Proton emission in both spherical and deformed systems has been studied extensively 
in the past decades Qj. For the spherical proton emitter, a simple WKB estimation of the 
transmission through the Coulomb and centrifugal barriers could give the correct order of 
magnitude of the decay rates and the angular m0m entn m of the decaying state g Q,B,Q. 
Most of the proton emitters in the rare earth region, which are predicted to have large static 



quadrupole deformations [3] are analyzed by a particle-coupled core model with the unbound 
proton interacting with an axially symmetric deformed core js, 9, 



id, 



11 



12 



3,Q, 



15 



16 



171 ]. Such an analysis over the past several years turned out to be a good description of the 



ground-state properties of axially deformed rare-earth proton emitters. 

. lecently, the proton emission from triaxial nuclei has drawn lots of attention 
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2d, 



221 ] . Specific combinations of single-particle orbitals near the Fermi surface can lead to 



a propensity toward triaxial shapes as illustrated by recent calculations of the additional 



binding energy due to non-axial degrees of freedom f 
of triaxiality throughout the nuclear chart. In Ref. Ij]] 



231 ] . which revealed several "islands" 
the static triaxial deformation was 



introduced in the adiabatic coupled channels method in order to investigate proton emission 
from 141 Ho(7/2~). The total decay rate and the 2+ branching ratio, however, were found 
to be good agreement with experimental data only for the triaxial angle 7 < 5°. The 
importance of triaxial deformation in proton emitters 161 Re and 185 Bi was pointed out in 
Ref. 20|], where the decay widths were found to be very sensitive to the 7 deformation. 
However, the sensitivity of decay widths in 161 Re and 185 Bi on the triaxial deformation 
was questioned in a recent paper reporting a non-adiabatic quasipartilce calculation 21]. 
Instead, the pairing effect was found to have a more significant influence. In Ref. 22J, the 
quasiparticle-coupled core model has been used to address the important role of the gamma 



degree of freedom in the prediction of the proton decay rate and the spectrum of excited 



3 



states in the proton emitter 145 Tm. 

It should be pointed out that the nuclear potentials in the (cmasi)particle-coupled core 
model are tuned to fit the measured energy of the decaying state 16(. In particular, it cannot 
provide any information about the microscopic structure properties of proton emitters. 

Approaches based on concepts of non-renormalizable effective relativistic field theories 
and density functional theory provide a very reliable theoretical framework for studies of 
nuclear structure phenomena at and far from the valley of /^-stability. In particular, the 
relativistic mean field (RMF) theory, which can take into account the spin-orbit coupling 
naturally, has been successfully applied in analysis of nuclear structure over the whole peri- 



odic table, from light to superheavy nuclei with a few universal parameters [24!.l25l.l26l.l27l.l28|. 
The application of RMF theory with the restriction of axial symmetry to study the proper- 
ties of proton emitter has already been done 



29 



30 



31] . In general, the predicted location 



of the proton drip line, the ground-state quadrupole deformations, one proton separation 
energies at and beyond the drip line, the deformed single-particle orbit occupied by the odd 
valence proton, and the corresponding spectroscopic factor are in good agreement with the 
experimental data. However, the influence of triaxiality on proton emitters has not been 
investigated in the microscopic self-consistent RMF approach. Here in this paper, the influ- 
ences of 7 deformation degree of freedom and pairing correlations on proton emitters 145 Tm 
will be studied within the triaxial deformed RMF approach. 

The paper is arranged as follows. In Sec. HIl a brief introduction of the RMF approach 
will be given. Both axial and triaxial calculations with pairing correlation are carried out 
to investigate the properties of proton emitter 145 Tm in Sec. IHIl including total energy, 
quadrupole deformations, one-proton separation energy, the potential of the valence proton, 
and the corresponding spectroscopic factor. Finally, our conclusions and summary are given 
in Sec. ED 



II. THE RELATIVISTIC MEAN FIELD THEORY 

The starting point of the RMF theory with meson-exchange providing nucleon-nucleon 
interaction is the standard effective Lagrangian density constructed with the degrees of 
freedom associated with the nucleon field (?/>), two isoscalar meson fields (a and u^), the 
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isovector meson field and the photon field (A^), 



i^d p -m- g a o - g^uj^ - g p ^r ■ p p - -e(l - t^A^ 



l -d»ad,o - U a (a) - ifi^fi^ + Ef w (w M ) - ^ • ^ + t^) 
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(1) 



where m and rrii(gi) (i = a, p M ) are the masses (coupling constants) of the nucleon and 
the mesons respectively and 



F" u = d»A v - d v A tx 



(2a) 
(2b) 
(2c) 



are the field tensors of the vector mesons and the electromagnetic field. Here in this paper, 
we adopt the arrows to indicate vectors in isospin space and bold types for the space vectors. 
Greek indices /i and v run over 0, 1, 2, 3, while Roman indices i, j, etc. denote the spatial 
components. 

The nonlinear self-coupling terms U a (cr), U^uj^), and U p (p^) for the a-meson, cu-meson, 
and p-meson in the Lagrangian density ([1]) respectively have the following forms: 



Ua(a) = \rnla 2 + ^g 2 a 3 + -g 3 a 4 , 
Uuiup) = -mf^o;,, + -c 3 (o/ 1 ^) 2 



-m 2 p p M ■ fa. 



(3a) 
(3b) 
(3c) 



In the mean-field approximation, the correspondent energy functional is obtained as 



^RMF \Pi 



d xTr 



(3[~f-p + m + g a a + g^uf + grfr ■ p p + -e(l - r 3 ) 7 ^ t ) p 



+ J d 3 x S^d°ad a - ^Padta + U a (a) - ^Q ^ + ~ U„ 



(4) 



where m denotes {a,u p , p^, A^} respectively. 
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The equations of motion for the nucleon and the mesons can be obtained by requiring 
that the energy functional (jlj) be stationary with respect to the variations of p and <f> m , 



5 {E RMF [p, 



Tr(ep)}=0, 



(5) 



where e is a diagonal matrix, whose diagonal elements are the single particle energies. Using 
the variation Sp with respect to ipk, the stationary condition (jSJ) leads to the Dirac equation, 



[a. ■ p + p(m + S + 'fVjtyk = e k ip k 



(6) 



for the nucleon. The scalar potential S and vector potential U M in Eq. (jfJJ) are respectively, 



S = g<rcr, 

Vf, = g^f, + g p r ■ p p + -e(l - r 3 ) A M . 
The Klein-Gordon equations for the mesons and the photon are given by, 



0, 



+ U'( 



±s, 



(7a) 
(7b) 

(8) 



where the (+) sign is for vector fields and the (— ) sign for the scalar field. The source terms 
£^ m in Eq.® are sums of bilinear products of Dirac spinors 



fc>0 
k>0 
k>0 



^ fc>0 



m = 0" 

m P/j, 



(9) 



where the sums run over only the positive-energy states (k > 0) (i.e., no sea approximation) 
and the occupation probability of the single-particle energy level k, i.e., v%, is evaluated 
within the BCS method. It is sufficient for proton-rich nuclei because of the attenuating 
effect of the Coulomb barrier on the spatial extent of the proton wave function, which limits 
the formation of proton halo. 



III. RESULTS AND DISCUSSION 

The static Dirac equation (jfJJ) for the nucleon and Klein-Gordon equation (jSj) for the meson 
fields are solved by expansion on the cylindric (axially deformed RMF) or three-dimensional 
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Cartesian (triaxially deformed RMF) harmonic oscillator basis with major shell numbers as 
Uf and nb for the nucleons and mesons respectively. The equation of motion for the photon 
field is solved using the standard Green's function method because of its long range. The 
parameter set PK1 33] is used throughout the calculation and the center-of-mass (cm.) 
correction is taken into account by 

= ds^- (10) 

where P c . m . is the total momentum of a nucleus with A nucleons. In order to check the 
convergence of the results with the number of expanded oscillator shells for nucleons rif and 
for mesons rib, the total energy, quadrupole deformation (3 and 7 of 144 Er as functions of shell 
numbers are calculated with axially (left panel) and triaxially (right panel) deformed RMF 
approaches as shown in Fig. [TJ It indicates that as long as rif > 14, rib > 18, the binding 
energies and the deformations are independent of the expanded shell numbers. Therefore 
in the following, rif = 14 and rib = 20 will be adopted. More details about the axially and 
triaxially deformed RMF approaches can be found in Ref. 34j and Refs. 35|, |36j respectively. 



In order to get the pairing gaps for Er and Tm, we fit odd-even mass differences of 



around 120 nuclides 



38 



in the very proton-rich side ranging from La to Re by the four-point 
difference formula. The neutron and proton pairing gaps obtained are, A^ st ' = 13.7/\/]4 
and Ap Sti = 15.9/v^4 respectively, and the corresponding rms deviation with respect to 
experimentally known values are 0.17 and 0.16 MeV. Considering the blocking effect of odd 
valence proton, the proton pairing gap for 145 Tm is reduced by a factor /, ranging from zero 
to one in the calculation. For the odd-mass system, the blocking calculations are performed 
without breaking the time-reversal invariance. In this case, the space-like components of the 



36|, 



vector fields vanish and this may introduce uncertainty around several hundred keV 
which may be compensated by the pairing gap. 

The total energy, deformation parameters f3 and 7 in 145 Tm as functions of the pairing 
reduction factor / are investigated in the triaxial RMF+BCS/PK1 approach and plotted 
in Fig. |2l in which the total experimental energy for 145 Tm is obtained from the systematic 
estimated atomic binding energy in Ref. 38[ by subtracting the electron binding energy 
according to Ref. (39]. 

It shows that the proton pairing gap A p has a significant influence on the total energy, 
but a negligible influence on the predicted deformation parameters. Especially the triaxility 
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parameter 7 is almost independent on A p . The one proton separation energy S p in 145 Tm as a 
function of the / is plotted in Fig. [31 It is noticed that the S p in the triaxial RMF+BCS/PK1 
calculation coincides with the experimental value, 1.728(10) MeV [3?! with a certain A p 
(i.e., / ~ 0.9). However, the axial RMF+BCS/PK1 calculation can not reproduce the 
experimental data of S p . It indicates the importance of both the 7 degree of deformation 
and pairing correlations in the description of proton emission in 145 Tm. 

The one proton separation energy S p , charge radius r c , neutron radius r n , as well as 
deformation parameters (3 and 7, the single-particle orbital occupied by the odd valence 
proton, and the corresponding spectroscopic factor u\ for 145 Tm are calculated in the ax- 
ially deformed and triaxial RMF+BCS/PK1 approaches ( / = 0.9 for proton in 145 Tm). 
The spectroscopic factor Sk of the deformed odd-proton orbital k is approximately given by 
the unoccupied probability u\ of state k in the daughter nucleus with an even proton num- 
ber 30] . The results are compared with those of relativistic Hartree-Bogoliubov (RHB), the 



Hartree-Fock-Bogoliubov (HFB-14), the finite-range droplet mass model (FRDM) as well as 
the experimental data in Table [H Similar to the HFB-14 prediction, the axially deformed 
RMF approach predicts an oblate shape with /3 = —0.21 for 145 Tm, while the RHB shows a 
prolate ground-state. The different predicted shapes for 145 Tm can be ascribed to the shape 
coexistence 3l|. While the large difference for one proton separation energy between the ax- 



ially deformed RMF+BCS/PK1(S , P = 0.62 MeV) and the RHB calculation (S p = 1.43 MeV) 
is ascribed to the treatment of the pairing. 

After taking into account the 7 degree of deformation self-consistently, the one proton 
separation energy (S p = 1.71 MeV) in the triaxial RMF+BCS/PK1 calculation is in good 
agreement with the data. The corresponding deformation parameters are respectively f3 = 
0.22 and 7 = 28.98° as shown in Table [H Furthermore, the spectroscopic factor of the odd 
valence proton is 0.67, which is consistent with the value 0.51(16) obtained with the WKB 

n 

approximation calculation [371 ] . 

The triaxiality can make the proton tunneling through the Coulomb barrier easier. The 
mean-field potential and density distribution of protons in 145 Tm are plotted in Fig. [4] 
as functions of x (for y = 0.52 fm and z = 0.52 fm) (dotted line), y (for x = 0.52 fm 
and z = 0.52 fm) (dashed line), as well as z (for x = 0.52 fm and y = 0.52 fm) (solid 
line), respectively. It shows that both the potential and density distribution are triaxially 
deformed. The Coulomb barrier is different in different directions and obviously the Coulomb 
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barrier in the direction y is lower than those in x and z. 

Apart from the Coulomb barrier, the proton decay probability also depends on its energy 
and orbital angular momentum, i.e., the centrifugal barrier. In Fig. [5J the single-particle 
energy levels for both neutron and proton in 145 Tm are given, in which each level is labeled 
with the quantum numbers of its main-component in the spherical Dirac spinor. The valence 
proton belongs to the /in/2 subshell, which is consistent with the observed spin-parity in 
Ref. 

In refs. 



41 



421 ] . fine structure in proton emission was observed in 145 Tm. In order to 
reproduce the experimental partial proton half-lives, it is found that the wave function of 
the valence proton in 145 Tm is composed mainly of 67% for 0/in/ 2 and 3.7% for I/7/2, which 
coupled to the ground state and the excited state of the 144 Er core. Using the particle-core 
vibration model 43j and assuming deformation (3 = 0.18, the experimental half-life and the 



fine structure branching ratio can be reproduced, in which the wave function of the valence 



proton is composed of 56% for O/in/2 and 3% for I/7/2 [42]. In contrast, by reproducing the 



fine structure branching, the particle-core vibration calculations in ref. 44j give 33% only 



for O/in/2- 

It is interesting to examine the composition of the valence proton in 145 Tm obtained from 
the present microscopic and self-consistent RMF calculation. In Fig. [6], the main spherical 
components of wave function for the valence proton in 145 Tm by axially and triaxially de- 
formed RMF calculations are given. The main components in axial RMF calculation are 
88.6% for O/in/2, 3.7% for O/19/2 and 2.3% for lf T/2 . While in triaxial RMF calculation, the 
main components are 82.1% for Ohu/2, 7.9% for I/7/2, and 2.3% for O/7/2. 

The potential and density distribution for the valence proton are obtained in triaxial 
RMF calculation. To show what the triaxial potentials and proton distribution look like, 
we plot the potential and density distribution for the valence proton in Fig. [71 in which the 
potential Vk(r) is given by, 



V k (r) = V (r) + S(t) + J^E WW + £ |GV| a V,(r)], Vi(r) = ^J^, (11) 



0„K. 



where and are respectively the expansion coefficients of the large and small com- 
ponents in spherical basis = \n£jrrij). 

A triaxial RMF+BCS/PK1 calculation with the estimated pairing gaps is also carried 
out for the neighboring proton emitters 146 Tm and 147 Tm. In contrast with the predicted 



shape transition from prolate to oblate existing from 145 Tm to 146 Tm [t[ H], large triaxial 
deformations have been found for 146 Tm (/3 = 0.19, 7 = 39.70°) and 147 Tm (/3 = 0.21, 
7 = 28.16°). The calculated one proton separation energies are 1.26 MeV and 0.90 MeV, 
which is close to the data, 1.120(10) MeV [45| and 1.054(19) MeV [46| respectively. 



IV. SUMMARY 



The axially and triaxially deformed RMF approaches have been applied for the description 
of ground-state properties of the proton emitter 145 Tm with parameter sets PK1 and a 
constant pairing gap BCS-method for the pairing correlations. It has been found that the 
triaxiality and pairing correlations are essential to reproduce the one proton separation 
energy in 145 Tm. The observed spin and parity of the emitted proton can be understood 
from the main spherical component by the transformation of the quantum number of the 
valence proton. The corresponding spectroscopic factor in 145 Tm obtained in the present 
calculation is consistent with that obtained with the WKB approximation calculation. Large 
triaxial deformations have been found in 146 Tm and 147 Tm as well. 
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TABLE I: Ground-state properties of the proton emitter Tm. The one proton separation energy 
S p , charge radius r c = \Jrp + 0.64, neutron rms radius r n , deformation parameters (3 and 7, as well 
as the quantum numbers of its main-component in the spherical Dirac spinor j n and spectroscopic 
factor u\ for the valence proton in triaxial RMF calculation with PK1, and pairing correlations 
treated by the BCS approximation are compared with those of the axial RMF calculation with 
PK1, and pairing correlations treated by the BCS approximation (Axial), the relativistic Hartree- 
Bogoliubov calculation (RHB), the Hartree-Fock-Bogoliubov calculation (HFB-14), the finite-range 
droplet mass model (FRDM) and experimental data. The neutron and proton pairing gaps in 144 Er 
and 145 Tm are estimated by A^ st - = 13.7/-/A MeV, A c p st - = 15.9/-/A MeV while the proton pairing 
gap in 145 Tm is given by 0.9 x 15.9/v^4 MeV. 
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FIG. 1: (Color online) The total energy, deformation parameters (3 and 7 calculated in axially (left 
panel) and triaxially (right panel) deformed RMF with PK1 for 144 Er as functions of the number 
of expanded oscillator shells for the meson field and the nucleon field rif. 
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FIG. 2: The total energy, deformation parameters (3 and 7 as functions of the proton pairing 
reduction factor in 145 Tm calculated by the triaxial RMF approach with PK1, where a reduction 
factor / = A p / Ap St ' varying from to 1 has been used in the BCS approximation and the data is 
given by a grey line. 
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FIG. 3: One proton separation energy plotted as a function of proton pairing reduction factor 
varying from to 1 for 145 Tm in the axial and triaxial RMF calculations with PK1, and pairing 
correlations treated by the BCS approximation. The neutron pairing gaps for 145 Tm and 144 Er as 
well as the proton pairing gap for 144 Er are chosen as the estimated values. The data is given by 
a grey line. 




FIG. 4: The mean-field potential and density distribution for the protons plotted as functions of 
x (for y = 0.52 fm and z = 0.52 fm) (dotted line), y (for x = 0.52 fm and z = 0.52 fm) (dashed 
line), as well as z (for x = 0.52 fm and y = 0.52 fm) (solid line) in 5 Tm. 
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FIG. 5: (Color online) Single-particle levels for neutrons and protons in 145 Tm in the triaxial RMF 
calculation with PK1, and pairing correlations treated by the BCS approximation. The level is 
labeled with the quantum numbers of its main-component in the spherical Dirac spinor. The black 
(blue) one is the level with positive (negative) parity. 
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FIG. 6: (Color online) Composition of the wave function for the valence proton in 145 Tm calculated 
with axially and triaxially deformed RMF approaches. 




FIG. 7: The proton-nucleus potential (upper panel) and density distribution of the valence proton 
(lower panel) plotted as functions of x (for y = 0.52 fm and z = 0.52 fm) (dotted line), y (for 
x = 0.52 fm and z = 0.52 fm) (dashed line), as well as z (for x = 0.52 fm and y = 0.52 fm) (solid 
line) in triaxial RMF calculation with PK1 and pairing correlations treated by BCS approximation. 
The proton Fermi level is given by the short-dotted line (upper panel). 
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